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"^t ! ABSTRACT 

Context. Pulsars PSR J0248+6021 (with a rotation period P = 217 ms and spin-down power E = 2.13 x 10 35 erg 
s" 1 ) and PSR J2240+5832 (P = 140 ms, E = 2.12 x 10 35 erg s" 1 ) were discovered in 1997 with the Nangay radio 
telescope during a northern Galactic plane survey, using the Navy-Berkeley Pulsar Processor (NBPP) filter bank. The 
GeV gamma-ray pulsations from both were discovered using the Fermi Large Area Telescope. 

Aims. We characterize the neutron star emission using radio and gamma-ray observations, and explore the rich envi- 
ronment of PSR J0248+6021. 

Methods. Twelve years of radio timing data, including glitches, with steadily improved instrumentation, such as the 
Berkeley-Orleans-Nangay pulsar backend, and a gamma-ray data set 2.6 times larger than previously published allow 
detailed investigations of these pulsars. Radio polarization data allow comparison with the geometry inferred from 
gamma-ray emission models. 

Results. The two pulsars resemble each other in both radio and gamma-ray data. Both are rare in having a single 
gamma-ray pulse offset far from the radio peak. The anomalously high dispersion measure for PSR J0248+6021 (DM 
= 370 pc cm -3 ) is most likely due to its being within the dense, giant HII region W5 in the Perseus arm at a distance of 
2 kpc, as opposed to being beyond the edge of the Galaxy as obtained from models of average electron distributions. Its 
large transverse velocity and the low magnetic field along the line-of-sight favor this small distance. Neither gamma-ray, 
X-ray, nor optical data yield evidence of a pulsar wind nebula surrounding PSR J0248+6021. We report the discovery 
of gamma-ray pulsations from PSR J2240+5832. We argue that it could be in the outer arm, although slightly nearer 
than its DM-derived distance, but that it may be in the Perseus arm at half the distance. 

Conclusions. The energy flux and distance yield a gamma-ray luminosity for PSR J0248+6021 of L 1 — (1.4±0.3) x 10 34 
erg s" 1 . For PSR J2240+5832, we find either L 7 = (7.9 ± 5.2) x 10 34 erg s^ 1 if the pulsar is in the outer arm, 

or L y = (2.2 ± 1.7) x 10 34 erg s _1 for the Perseus arm. These luminosities are consistent with an L 7 cx \[~E rule. 
Comparison of the gamma-ray pulse profiles with model predictions, including the constraints obtained from radio 
polarization data, implies outer magnetosphere emission. These two pulsars differ mainly by their inclination angles 
and acceleration gap widths, which in turn explain the observed differences in the gamma-ray peak widths. 
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1. Introduction 

Pulsars PSR J0248+6021 and PSR J2240+5832 were dis- 
covered with the Nancay radio telescope (NRT) in a sur- 
vey of the northe r n Galactic plane from 1997 to 1999 
(jFoster et al.lll997ft . iRav et all (|1999D reported their loca- 
tions and rotation periods. Twelve years of timing measure- 
ments allow us to provide accurate ephemerides, including 
a proper motion determination for PSR J0248+6021 that 
helps us understand its relation to the complex of gas clouds 
co-located on the sky. Polarization data provides geometry 
constraints in addition to those obtained by comparing the 
gamma-ray pulse profiles with model predictions. 

By coincidence, both pulsars have a high spin-down 
power E = 2.1 x 10 35 erg s -1 , making it als o likely that 
they e mit in GeV gamma-rays, as discussed bv lSmith et al.l 
(|2008h . Indeed, PSR J0248+6021 is among the 46 gamma- 
ray pulsars descr ibed in the "First Pulsar Catalog" 
(jAbdo et al.ll2010bl) using the Large Area Telescope (LAT) 
on the Fermi Gamma-ray Space Telescope (formerly 
GLAST). They have the same magnetic field strengths at 
the neutron-star light cylinder, Blc, to within 20%, and 
the characteristic ages t c = P/2P differ by only a factor 
of two. Table [1] lists measured and derived parameters for 
both pulsars. PSR J2240+5832 has a smaller dispersion 
measure (DM = 263.5 pc cm~ 3 ) than PSR J0248+6021 
(DM = 370 pc cm- 3 ), and t he NE2001 model o f the 
Galactic electron distribution (|Cordes k Laziol |2002[) as- 
signs PSR J2240+5832 a distance of 10.3l~ 3 kpc, whereas 
it places PSR J0248+6021 well beyond the edge of the 
Galaxy (> 43.5 kpc). Curiously however, it is the nominally 
closer pulsar, PSR J2240+5832 that pulsates less brightly 
in the Fermi data. 

Two obvious solutions to this seeming paradox are that 
either there is some intrinsic difference between the pul- 
sars, such as their orientation angles and beaming, or that 
the distances are simply misunderstood. In this paper, we 
clarify the situation. 



2. Nancay northern Galactic plane survey 

The NRT is a meridian instrument, with a primary an- 
tenna 200 m wide by 35 m tall, equivalent to a parabolic 
dish with a 94 m diameter. It can track objects with de- 
clinations S > —39° for roughly one hour around culmina- 
tion. The half-power beam width at 1.4 GHz is 4 arcmin 
(east-west) by 22 arcmin (north-south) at S = 0°, a shape 
well-adapted to sky scanning. The 1.4 GHz nominal system 
temperature at the time of these observations was about 50 
K in both horizontal and vertical polarizations, for S ~ 60°. 
The "FORT" receiver upgrade from 1999 to 2001 improved 
sensitivity by a factor of 2.2, with an efficiency of 1.4 K/Jy 
and a system temperature of 35 K at 1.4 GHz. The fre- 
quency coverage is now continuous from 1.1 GHz to 3.5 
GHz and allows full Stokes measurements (jTheureau et al.l 
l2005h . 

Two pulsar instrumentation systems ( "backends" ) were 
used successively, first for the survey and then for the 
continued timing. The earlier system was the Navy- 
Bcrkclcy Pulsar Processor (NBPP) filter bank, designed 
and built for this survey at the Naval Research Laboratory 
(NRL) in collaboration with the University of California, 
Berkeley. We set the 96 channels in the full crate to 
cover 1.5 MHz/channel, for a total bandwidth of 144 MHz 



([Backer et a,l.lll997h . The output data has a time resolution 
of 50 to 100 /is in search mode, allowing detection of signals 
with periods as short as 0.1 ms. For the survey, data were 
acquired at a sampling rate of 60 /zs using 4-bit digitization 
of 96 channels and both polarizations. 

Since 2004, the NRT pulsar backend is the BON 
(Berkeley-Orleans-Nangay) coherent dedispersor with a 
bandwidth of 128 MHz. A spectrometer digitizes data 
voltages, followed by four data servers that share the 
data over a 70-node cluster of personal computers (PCs) 
running the Linux operating system. In 2009, a GPU- 
based (Graphics Processor Unit) computer system replaced 
the PC cluster, with the same computing power in only 
two nodes. Dedispersion is performed with a special fil- 
ter in the complex Fourier domain. The timing resolution 
i s a few 100 ns on the most stable millisecon d pulsars 
(jCoenard k Theureaull200l ICoenard et al.ll2009h . 
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Fig. 1. Phase-aligned 7-ray and radio light curves for PSR 
J0248+6021 obtained with the Fermi Large Area Telescope 
and the Nangay Radio Telescope. The bottom panels show 
the radio profiles at three frequencies used to build the 
ephemeris. The second panel from the bottom shows the de- 
gree of linear (red dashed) and circular polarizations (blue 
dotted), as well as the linear polarization position angle and 
a Rotating Vector Model (RVM) fit. The other panels show 
the phase- folded 7-ray data in different energy bands. Two 
rotations are shown for clarity. 
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Fig. 2. Phase-aligned light curves for PSR J2240+5832. 
The bottom panel shows the 1.4 GHz radio profile, along 
with the linear (red dashed) and circular polarizations (blue 
dotted). The second panel from the bottom shows the lin- 
ear polarization position angle sweep with an RVM fit. The 
other panels show the phase-folded 7-ray data obtained 
with the Fermi Large Area Telescope in different energy 
bands. Two rotations are shown. 



The NBPP was used to conduct a high-frequency, wide- 
bandwidth (1280-1430 MH z) pulsar survey of the Galactic 
plane from 1996 to 1998 ([Foster et all 1 1996(1 . The search 
covered longitudes —15° < I < 180° for Galactic latitudes 
|6| < 3°, for a total sky coverage of 1170 square degrees. 
The survey was optimized to find distant pulsars with peri- 
ods as short as the theoretical break-up speed of a neutron 
star. Over 40 000 pointings were acquired, each with an 
integration time of 2 minutes, for an expected 5tr detec- 
tion of a 0.5 mjy peak intensity profile after dedispersion. 
Observations were completed in late 1998. The final data 
volume is ~ 5 Tbytes. 



A total of two pulsars were discovered in the survey, during 
the first 1998 data analysis based on a fast Fourier trans- 
form (FFT). Figs. [T]and [2] show multiwavelength pulse pro- 
files for PSRs J0248+6021 and J2240+5832, respectively. 
Table[T]lists various properties. PSR J0248+6021 has a spin 
period of 217 ms, with a duty cycle of 12% at 1.4 GHz. A 




Fig. 3. Expanded view of the radio polarization position 
angle sweep near the peak in radio intensity. The red points 
show the data used in the RVM fit. The black points failed 
the selection cuts described in the text. Top two frames: 
For PSR J2240+5832, at 1.4 GHz. The RVM curve shown 
corresponds to inclination angles a = 108° and ( = 123°. 
Bottom two frames: For PSR J0248+6021, at 2.1 GHz. The 
RVM curve shown corresponds to inclination angles a = 
46° and C = 52°. 
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Fig. 4. A single pulse search of PSR J0248+6021. For each 
of 250 DM values, the 96 NBPP radio frequency samples are 
de-dispersed and summed. Bottom frame: Summed inten- 
sities with signal-to-noise 5<r above the average noise level 
are shown as points. Larger points indicate larger S/N. Top 
frames: Projections showing a clear pulsar detection. 



4 



Theureau et al.: PSRs J0248+6021 and J2240+5832 



100 



PSR J0248 + 6021 «=-1.19±0.07 



-,10 



PSR J2240 + 5832 a = - 1 .83±0.25 



0.1 



Frequency (GHz) 

Fig. 5. Calibrated mean radio flux densities as a function 
of observation frequency. Note that the sign convention for 
the index is the opposite of that used for the gamma-ray 
spectra. 



scattering tail is seen at low frequency, as expected from 
the large DM and confirmed by the observations made at 
higher frequencies. PSR J2240+5832 has a spin period of 
140 ms and the radio peak is even narrower, 7% of a ro- 
tation at half-maximum. Relatively few pulsars are known 
in this direction (/ = 106.6°) and none with such a large 
DM. It may be one of the few pulsars known in the outer 
arm of the Galaxy and amongst the most distant known 
gamma-ray pulsars. 

Subseque nt FFT re-analysis us ing the PRESTO package 
developed bv lRansom et al.l (|2003|) identified no additional 
pulsars. We also used PRESTO to search for single pulses 
from the entire dataset: raw data were first dedispersed for 
250 values ranging from to 500 pc cm -3 keeping the origi- 
nal 60 /zs resolution. Each time sequence was then searched 
for single pulses with a signal-to-noise ratio (hereafter S /N) 
greater than 5. The search used matched filtering with box- 
cars up to 30 bins allowing good sensitivity up to 1.8 ms 
pulse width. We clearly detected series of very strong bursts 
from PSR J0248+6021 at the pulsar DM as shown in Fig.ffl 
No other single pulse signal was seen, including from PSR 
J2240+5832. PSR J0248+6021 is thus somewhat rare, in 
that it seems to burst in a way similar t o PSR 130656+14 
and R RATs (Rotating Radio Transients, IWeltevrede et all 
l2006t ). The single pulse search is confirmed to be comple- 
mentary to a classical FFT analysis, favoring pulsars with 
a low average signal strength but having strong bursts. 



3. Radio measurements 

3.1. Flux density 

The radio flux for pulsar profiles obtained with the Nangay 
coherent dedispersion instrumentation can be accurately 



measured using a pulsed noise diode. The diode's counts-to- 
mjy scaling is obtained from fiducial galaxies. The diode 
is fired at 3.3 Hz for 15 seconds at the start of each ob- 
servation, then phase-folded as for a pulsar. Several dozen 
observations at 1.4, 2.05, and 2.68 GHz were selected for 
PSR J0248+6021. The 1.4 GHz average value for the mean 
(i.e., phase-averaged) flux density of 13.7 ± 2.7 mJy is 
shown in Fig. [51 and a spectral index of —1.19 ± 0.07 is 
obtained. These values are consistent with those obtained 
previously using the radiometer equation for both polariza- 
tion channels, along with the known telescope characteris- 
tics (T sys = 20.5 Jy) and an assum ed continuum flux of 6.5 
Jy at 1.4 GHz in this sky direction (Reich 1982), to convert 
rms noise fluctuations in the off-pulse part of the profile to 
a Jansky scale. 

Fourteen observations for PSR J2240+5832 yield a 1.4 
GHz average value for the mean flux density of 2.7 ± 0.7 
mJy, and a spectral index of —1.83 ± 0.25, as also shown 
in Fig. [3] The sign convention for the radio energy spectral 
index is the opposite of that used for the gamma-ray photon 
index (i.e., S(v) oc v a as compared to E~ r as in Equation 
2 in Section 5). 



3.2. Polarization 

Fig. [3] zooms in on the sweep with phase of the linear po- 
larization position angle (PA) for both pulsars. PA data 
is shown for all points with a minimum S/ N in Figs. [1] 
andEl We reduced the data using PSRCHIVE (jHotan et all 
2004). For PSR J0248+6021, we averaged 399 minutes of 
data from MJD 54851 to 55182, and obtained a rotation 
measure of RM = -158 ± 6 rad m~ 2 (3a statistical uncer- 
tainties). At inflection, extrapolated to infinite frequency, 
we obtain PAo = —67° ± 3°, shown in Fig. [6] together 
with the proper motion direction, discussed below. For PSR 
J2240+5832, we used 846 minutes of 1.4 GHz data recorded 
from MJD 54866 to 55209, yielding RM = 24 ± 4 rad m~ 2 
and PA = -39° ± 2°. 

We fit the PA and phase data with the ro tating vector 
model (RVM, Radhak rishnan fc Cookd Il969h to estimate 
the inclination £ of the neutron-star rotation axis to the 
line-of-sight to Earth, and the angle a between the rotation 
axis and the magnetic axis. We use only data points with 
a minimum S/N of cr PA < 2°. For PSR J0248+6021, the 
PA sweep inflection point seems to lead the radio peak due 
to bias induced by the broad scattering tail prominent at 
1.4 GHz. We therefore use the 2.1 GHz data for the RVM 
analysis. The fit results still vary with the choice of data in 
the tail but become stable when we favor the leading edge 
of the PA swing. We interpret this as a deformation due 
to residual scattering. We reject the remaining points with 
phase 0.01 beyond the peak. The results are robust once 
the choice to cut the tail is made. The points used for the 
RVM fit are shown in red in Fig. [3] As usual for RVM fits, 
the magnetic impact parameter f3 — £ — a is most strongly 
constrained, here to j3 ~ +5°. A wide range of a from 40° 
to 80° provides comparably good fits, and acceptable fits 
extend from a w 25° to 110°. The best fit has \ 2 = 86.6 
(reduced x 2 = 1-6). 

The narrow radio pulse of PSR J2240+5832 is less af- 
fected by scattering: a fit to the 1.4 GHz profile yields 
/? = 16° (for typical a). Dispersion in the trailing data 
points causes a poor fit, with a minimum at x 2 = 329 (re- 
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Fig. 6. Environment of PSR J0248+6021 from red Palomar Optical Sky Survey photographic plates (POSS-II, back- 
ground) and NRAO VLA Sky Survey (NVSS) 1.4 GHz charts (contours). The solid arrows show the proper-motion 
position angle and uncertainty. The dashed arrows show the polarization position angle at the inflection point and un- 
certainty, extrapolated to infinite frequency. The arrow lengths have no meaning. The zoomed image at right shows the 
star near the pulsar excluded as an optical counterpart, as well as the optical arc emanating from the star discussed in 
the text. 



duced x 2 — 5.1). Comparable fits persist from a « 75° to 
130°, and plausible solutions extend from a ~ 10° to 150°.' 

The green contours in Figs. [7] and [8] summarize the 
RVM fit results, while the RVM curves in Fig. [3] show 
the fit for specific (a, £) combinations within the fa- 
vored region. These classic RVM fits assume radio emis- 
sion at low altitudes. It is becoming increasin gly clear 
(|Karastergiou fc Johnston] 120071 : iRavi et al.ll2010t ) that for 
many young pulsars, especially the gamma-ray detected 
pulsars, the radio emission arises at a substantial fraction 
of the light cylinder radius, Rlc- This in turn means that 
the polarizat ion PA sweeps are offset significantly from the 
radio pulse (Blaskicwi cz et al.lll991l ). and that more sub- 
tle distortions in the PA sweep shape may occur. Self- 
consistent fitting using large radio altitudes should there- 
fore lead to a shift in the inferred a and £ values. 



3.3. Radio timing and proper motions 





Date (MJD) 

Fig. 9. Time-of-arrival post-fit residuals for PSR 
J0248+6021 when glitch parameters are fit. Vertical 
lines indicate glitch epochs. The timing resolution of the 
datasets using the NBPP (left), BON prototype (middle), 
and BON (right) backends is apparent. 



Both pulsars were observed regularly for twelve years 
with increasingly advanced instrumentation, except for 18 
months in 2003-2004 when only the BON prototype with 
a 16 MHz bandwidth was available. Fig. EH shows the evo- 
lution of the timing precision over the years, and also in- 
dicates five principal glitch epochs for PSR J0248+2240. 
The largest was 5v/v = 7.5 x 10~ 7 in March 2009 (MJD 
54897.41), where v is the rotation frequency and 8v is 
the permanent pulse frequency increment. The epochs and 
8v/v values for the other four glitches were (52240.0, 1.6 x 
10- 10 ), (52850.0, 8.2 x 10" 9 ), (53700.0, 6.9 x 10~ n ) and 
(53900.0, -2.6 x 10" 11 ). Fig. [TU] shows the BON data, in- 
cluding the data used to phase-fold the gamma-ray pho- 
tons, in greater detail. The timing model includes a first 
time derivative of the DM ("DM1" variable in TEMP02). 
Higher-order DM variations as well as achromatic timing 



noise can be seen. For PSR J2240+5832, no large glitches 
were seen, but even after whitening some residual structure 
due to timing noise remains. 

Timing noise biases proper-motion measurement unless 
care is taken. Different approaches were compared, sum- 
marized in Fig. 1111 where the labels for the methods are 
described below. The first two methods compare the po- 
sition at different epochs with separate data spans, fit- 
ting for rotational parameters {v, z>, v) and pulsar posi- 
tion. The old NBPP timing data yield a position at MJD 
51250, with another position at MJD 53900 from the BON 
data, for a proper motion of \i a = 13.4 ± 8 mas yr _1 and 
Us = 58.2 ± 2.5 mas yr" 1 (label NBPP_2_BON). In a sim- 
ilar way, the BON data were divided into two equal parts. 
Two positions obtained at MJDs 53700 and 54500 indicate 
a proper motion of /j a = 42 ± 6 mas yr _1 and /15 = 55 ± 3 
mas yr _1 (splitted_BON). The two last methods use the 
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Fig. 7. Pulsar geometry and emission modeling for PSR J0248+6021. Green contours show the rotating vector model 
fit to the radio polarization data (contours at <5(x 2 /DoF) = +0.25, +0.5 above the minimum x 2 /DoF of 1.6). The color 
backgrounds are \3 maps of the fit to the observed > 100 MeV pulse profile to the outer gap model (left), the two-pole 
caustic model (middle), and the separatrix layer model (right), for different values of the magnetic inclination, a, and 
the minimum angle to the line-of-sight, ( (|Romani fc Wattersll2010l) . Each panel has the same color scale, where dark 
colors represent better fits. The preferred models lie along the green RVM-selected band. 




harmo nic whitening procedure proposed by iHobbs et al.l 
(|2004f) in their Appendix A. An iterative process using 
a fit of h armonically related sinusoids ("FITWAVE" in 
TEMP02 jEdwards et aT1l2006tl is applied following a strict 
sequence of different parameter adjustments. This proce- 
dure was applied to both the old NBPP and the new BON 
timing data. The NBPP data give fi a = 24 ± 30 mas yr" 1 
and us = 17 ± 12 mas yr" 1 (NBPP_FITWAVE) . The new 
BON timing data imply a proper motion of /i a = 62 + 3 
mas yr" 1 and fi s = 35.6+1.8 mas yr" 1 (BON_FITWAVE). 

The error bars shown in Fig. [IT] for the successive mea- 
surements sometimes do not overlap because systematic bi- 
ases due to e.g., timing noise or covariances underestimated 
by the fitting routines, are larger than the statistical uncer- 
tainties. We computed a weighted mean proper motion for 
PSR J0248+6021 based on all the determinations except 
the NBPP_FITWAVE, which has very large uncertainty 
and hence constrains little. The mean values fi a = 48 ± 10 
mas yr -1 and fis = 48 ± 4 mas yr -1 were obtained by 
weighting by the inverse of the uncertainties. The total is 
fi = \J p? a cos 2 (S) + = 53+11 mas yr -1 and the celestial 
position angle is PA = tan _1 (/^ Q cos(<5)/ fis) = 27° ± 6°. 



PSR J2240+5832's proper motion is roughly half that 
of the other pulsar. The measured value became sta- 
ble on ly after multiple iterations of the method cited 
above (Hobbs et al.l l2004f ). apparently due to one or two 
marginally detectable glitches that perturb the fits. We ob- 
tained /i Q = —6.1 ±0.8 mas yr -1 and fis = —21.0 ±0.4 mas 
yr -1 (^ = \J [i 2 a cos 2 (5) + /i 2 = 21.3 + 0.4 mas yr" 1 and ce- 
lestial position angle PA = 8.7° ±1°). The uncertainties are 
statistical. To be conservative, we assume systematic uncer- 
tainties of the same magnitude as those observed for PSR 
J0248+6021, and thus obtain fi — 21 ± 4 mas yr -1 (same 
fractional error) and celestial position angle PA = 9° ± 6° 
(same absolute value). 

Fig. [5] shows both the proper motion and polarization 
position angles for PSR J0248+6021 . They are rough l y per- 
pendicular (difference of 94° ±7°). Uohnston et~all |2007l) 
recall that neutron star velocity vectors and spin-axis may 
naturally become aligned at birth, if the supernova "kick" 
is preferentially along the spin-axis, and the radio emission 
is orthogonal to the magnetic field. That paper adds a new 
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Fig. 10. Time-of-arrival post-fit residuals for PSR 
J0248+6021, allowing for linear variations of disper- 
sion measure over time, as in Fig. |H1 but showing only data 
acquired with the BON backend. The 2.0 GHz data (blue 
squares) are clearly above the 1.4 GHz data (black dots) 
for the year 2006 and for 2009-2010, but below these data 
in 2008. We interpret the increasing, then decreasing DM 
value as evidence that the neutron star is moving through 
a dense environment. 
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Fig. 11. The PSR J0248+6021 proper motion measured 
by different methods (Section 3.3). The weighted average 
of /i = cos 2 (S) + JjR = 53 ± 11 mas yr _1 and PA = 
27° ± 6° is shown by the dark cross. The three different 
determinations are labeled NBPP_2_BON, splitted_BON, 
and BON_FITWAVE and are described in the text. Solid 
error bars are ±la and dashed error bars are ±3<r. 



sample of measurements for which this scenario appears to 
hold about half the time, weakening the observational evi- 
dence supporting this picture. For PSR J2240+5832, the 
difference between the velocity and polarization PA's is 
48° ± 6°: hence, in this new sample of two pulsars, again 
only half seem to comply with the "rule" . 



4. Distances and surroundings 

One goal of this work is to compare model predictions of 
the pulsars' gamma-ray luminosity L~, with o bservations. 



Follo wing the "Atlas" of IWatters et all (|2009h and its up- 
date (|Romani fc Wattersl 2010f ). we define 



4irfnG 100 d 2 



(1) 



gamma-ray flux radiated into all space to that along the 
line-of-sight, both averaged over a neutron star rotation. 
Most models depend on the inclination and viewing angles 
a, C as discussed below. The distance d is clearly important. 

4.1. Dispersion measure, kick velocities, and rotation 
measures 

4.1.1. Distance to PSR J0248+6021 

As stated earlier, the DM for PSR J0248+6021 is so large 
for this line-of-sight (I = 136.9°, b = 0.698°) that the 
NE2001 model places it outside the Galaxy. However, Fig. 
[5] shows extensive clouds that add to the electron column 
density and bias the NE2001 distance. 

The pulsar lies 0.3° west of the center of the open cluster 
IC 1848, which itself lies within the giant HII region W5, 
the "Soul" nebula in the Galaxy's Perseus ar m. The clus- 
ter diameter is 0.34° (jKharchenko et al.l[2005h . The pulsar 
was born [it c < 1° southwest of its current position. The 
heliocentric distance to IC 1848 estimated from both main 
sequence fitting of the star cluster and a Galactic kinemati c 
model converges towards 2 kpc (|Kharchenko et all [2005). 
That reference gives no distance uncertainty^]. Distance un- 
certainties for individual sta rs from other studies , e.g. using 
photometric measurements ([Morgan et al.ll2009l) do not ex- 
ceed ±0.4 kpc, and we assume ±0.2 kpc. We now argue that 
the pulsar is at the distance of the open cluster, i.e., on the 
near side of W5. 

The parallel component of the magnetic field along the 
line-of-sight to PSR J0248±6021 is B\\ = 1.23RM/DM = 
—0.5 /xG. Th i s mat ches nicely the expectation from Fig. 6 of 
iMitra et al.l ()2003[ ) and therefo re the pulsa r is m ost likely 
within the Perseus arm. From iReid et al.l (|2009), the far 
edge of that arm is < 3.6 kpc distant. 

The proper motion of PSR J0248+6021 is, with its un- 
certainty, > 42 mas yr . The largest measured neutron- 
star velocities due to su pernova "kicks" at birth are of the 
order of 1000 km s" 1 (|Hobbs et af]|2005l ). Thus, if PSR 
J0248+6021 were beyond 4.8 kpc it would have an unusu- 
ally large kick velocity. Conversely, a 2 kpc distance im- 
plies a very typical transverse velocity of vt ~ 500 km 
s -1 . Furthermore, the DM of PSR J0248±6021 increased 
steadily by 0.3 cm~ 3 pc per year from 2006 to 2008, and 
decreased at the same rate in 2009 and 2010. Fig. [TU] shows 
that the detailed behavior is complex. The simplest inter- 
pretation of the DM evolution is that the pulsar is moving 
rapidly in a dense environment, that is, that the pulsar is 
in a cloud (conceivably a nearer, invisible cloud). 

For the cloud coincident with the pulsar (Fig. [6]), the 
largest velocity (largest Doppler shift) of reliable HI fea- 
tures towards PSR J0248±6021 is -50 km s" 1 . In light 
of th e known velocity reversal in this part of the Perseus 
arm (jTian et al.ll2010f ) the cloud is thus between 2 and 3 
kpc, like the rest of W5. High resolution Ha imaging could 
confirm this preliminary indication. 

Towards PSR J0248+6021, the NE2001 model places 
the edge of the stellar Galactic disk 9 kpc from the Sun, 
with a DM of 210 cm -3 pc. A dense local environment, as 
within W5, may explain the large observed electron column 
density DM = 370 cm" 3 pc. For 2 kpc, NE2001 predicts 



for Gioo, the integral gamma-ray energy flux above 100 
MeV. The beam correction factor /q is the ratio of the 
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DM ~ 61 cm~ 3 pc, implying an excess of ~ 315 cm -3 pc. 
The main part of the HII region neighboring the pulsar is a 
shell of diameter 0.5° on the sky, equivalent to a true size of 
17.5 pc at 2 kpc. This converts into an electron density in 
the cloud of < n e >~ 18 cm~ 3 . Though large, this is typical 
of other HII regi ons, where values from . 1 to 1 00 cm ~ 3 
are seen (see e.g. iMitra fc Ramachandranl (|200lh . for their 
study of the scatter broadening in the direction o f the Gum 
Nebula). The WISE collaboration (|Wrightll2010l ) released a 
detailed image of the Heart and Soul nebulae seen between 
3.4 and 22 micronsfl PSR J0248+6021 lies in one of the 
"thickest" areas, where the shell is seen edge-on and the 
electron density is surely greater than the average value, 
indicating again that the large DM is compatible with a 
small pulsar distance. 

The large scattering tail excludes the notion that PSR 
J0248+6021 may be within the Orion spur, thus only ~ 1 
kpc away, because of nearby but invisible electron clouds. 
To conclude, PSR J0248+6021 is almost certainly in W5, 
at the distance of IC 1848, and we use d = 2.0 ± 0.2 kpc. 

4.1.2. Distance to PSR J2240+5832 

PSR J2240+5832 has DM = 263.5 cm" 3 pc and along 
this line-of-sight (Z = 106.57°, b = -0.111°) the NE2001 
distance is 10.3l^° 3 kpc. Archival optical and radio im- 
ages show no bright structures near the pulsar position 
that would indicate significant deviations from the NE2001 
maps. However, the pitch angle of the outer arm has been 
shown to be smaller than kn own when the NE2001 model 
was created (|Reid et al.ll2009D . and the far edge of the outer 
arm is < 8.3 kpc in this direction. 

Again assuming a transverse velocity vt < 1000 km s , 
the proper motion of > 17 mas yr" 1 of PSR J2240+5832 
limits its distance to < 12 kpc. Conversely, a value of 
vt = 400 km s _1 would place the pulsar in the Perseus 
arm, which extends from 3 to 4.6 kpc in this direction. PSR 
J2240+5832 is therefore among the few pulsars in this direc- 
tion for which rotation measures are available ([Han et al.l 
2006). The positive value of RM= 24 rad m" 2 , albeit small, 
is anomalous: other low-latitude RM's for 90° < I < 120° 
are negative. On the other hand, in this direction no other 
pulsar between 6 and 8 kpc has an RM measurement and 
the positive value may be due to the magnet ic field in- 
versio n beyond the Perseus arm suggested by lHan et ahl 
(j2006l) . The magnetic field B\\ = 1.23 RM/DM = 01 fiG is 
within the range of observed values in IMitra et~aT1 ((2003) . 

Neutron stars are most likely born within arms and, 
since vtt c < 0.1 kpc, the pulsar should not have gone far 
from its birthplace. We conclude by presenting two possi- 
ble values for the distance to PSR J2240+5832: one in the 
outer arm (d — 7.7 ± 0.7 kpc) and another in the Perseus 
arm (d = 3.8 ± 0.8 kpc). If in the outer arm, vt is large 
but not exceptional and the pulsar is farther than most 
known gamma-ray pulsars. If it is in the Perseus arm, the 
DM excess and positive RM remain unexplained but the 
speed and distance are more typical. We see below that the 
two values of L 7 obtained using the two distances are both 
plausible. VLBI astrometric parallax measurements would 
be of great interest for these and all gamma-ray pulsars. 



2 http:/ /antwrp. gsfc.nasa.gov/apod/apl00601. html 



4.2. Evidence of a pulsar wind nebula? 

PSR J0248+6021 coincides with the heart of an elongated 
nodule 15' long (se e Fig. H]), seen at both 21 cm (NVSS, 
ICondon et al.l [l998h and in the optical. The pulsar proper 
motion is ~ 45° from the nodule long axis. The POSS-II 
(Palomar Optical Sky Survey) image shows an arc 20" long 
about 5" from the pulsar position. The simple formulae 
for the size and offsets for pulsar wind nebula (PWN) ter- 
mination shocks given by iKargaltsev fe Pavlov! (|2008t ) are 
roughly consistent with what we observe. Might we be see- 
ing shocks in the interstellar medium caused by the pulsar's 
high transverse velocity? 

To explore t he idea of a po ssible PWN, we 
obtained Swift (jGehrels et al.l 120041) observations of 
PSR J0248+6021 with data taken over two days, 2009 Dec. 
10-11. The target was not d etected in the X-ray Telescope 
(XRT; iBurrows et aT1l2005D combined 7.2 ks exposure and 
we derive a 3<r limit on the source count rate (0.5-8 keV) 
of 0.0024 cts/s after dead time, vignetting, and PSF cor- 
rections. Using PIMMS, this corresponds to an observed 
flux limit of 1.4 x 10~ 13 erg cm~ 2 s _1 , or 2.0 x 10~ 13 erg 
cm" 2 s" 1 unabsorbed with Galactic Nh — 8.3 x 10 21 cm~ 2 
(iKalberla et al.l[2005h . assuming a power-law photon flux 
with T = 1.5 as is typical among observed X-ray PWN 
(same sign convention for V as in Eq. 2, Sect. 5). For d = 2 
kpc, we obtain (absorbed) Lx < 6 x 10 31 erg/s for the 
P WN X-ray luminosity Only e ight of the 40 X-ray PWN's 
in IKargaltsev fc Pavlovl (|2008h are dimmer: if there is an 
X-ray PWN, it is faint, possible for a PWN of the age of 
these pulsars. 

An upper limit on GeV PWN emission is given in 
Section 5.1. The position of PSR J0248+6021 is at the 
corner of b oth the MAGIC (lAlbert et al.l I2006D and the 
VERITAS (|Acciari et aH l2008ft TeV significance maps for 
LS 1+61° 303. No evidence of TeV emission from a putative 
PWN is apparent. Reanalysis of the complete TeV datasets 
optimized for a source at this location with PWN-like prop- 
erties (i.e. possibly offset, and/or spatially extended) is en- 
couraged. 

4.3. Upper limit for detection of an optical pulsar 

Digital Sky Survey images show a faint source seemingly 
coincident with the pulsar positio n. Analyzing the Ul tra- 
Violet/Optical Telescope fUVOT: Homing et all 12005] ) U- 
filter image (exposure: 3575 s) obtained with Swift on 2009 
December 11th, the faint source is detected at 6. Oct. We 
measure its center (J2000) at RA = 02h 48' 18.711", dec = 
60d 21' 38.98" with a statistical uncertainty of 0.46" (90% 
confidence). Its observed flux density is 1.07±0.18 xl0~ 2 
mJy at 8.56 x 10 14 Hz corresponding to magnitude 20.33 ± 
0.18. The uncertainties are statistical. This source is not 
detected in the UVOT Wl-filter image (3537 s) obtained 
on 2009 December 10th and we derive a 3cr upper limit 
of < 2.80 x 10~ 3 mJy at 1.14 x 10 15 Hz. Pulsar timing is 
incompatible with the position of the optical object and 
there is therefore no evidence of an optical pulsar. 

5. Gamma-ray Observations 

The LAT is an electron-positron pair conversion telescope 
that was placed in orbit with the Fermi satell ite on 2008 
June 11 (jAtwood et alJl2009l lAbdo et al.ll2009eD . The LAT 



Theureau et al.: PSRs J0248+6021 and J2240+5832 



9 



covers the 20 MeV to > 300 GeV energy range with 
higher sensitivity and more accurate localisation than pre- 
vious instruments (an on-axis effective area ~ 8000 cm 2 
above 1 GeV and angular resolution #68 ~ 0.8° at 1 GeV). 
LAT measurements yielded a catalo g of 1451 gamma-r ay 
"1FGL" sources during its first year (jAbdo et al.ll2010aT ). 

We used the standard Science Tools software package 
for Fermi LAT data analysifl For both pulsars, we se- 
lected data collected between 2008 August 4 (MJD 54682) 
when Fermi began scanning-mode operations, and the end 
of the ephemeris validity range (over 15 months in both 
cases). We kept "diffuse" class events (highest probability 
of being 7-ray photons) within a 15° "region-of-interest" 
(ROI) around the pulsar. We excluded events with zenith 
angles > 105° to reject the 7-ray albedo from the Earth's 
limb. Photon phases w ere calculated using the TEMP02 
pulsar timing software (Edward s et al.ll2006l) . 

At low photon energies, multiple scattering dominates 
the LAT's angular resolution. For the gamma-ray pulse pro- 
files (phase histograms), we applied an energy-dependent 
angular radius cut centered on the pulsar that approxi- 
mates the instrument point spread function (PSF), given 
by < e 6S (E) > = v /(5.12°) 2 (100 MeV/E) 1 - 6 + (0.07 ) 2 , 
which corresponds to a 68% containment angle (P6-V3 
Instrument Response Function, "IRF"). A maximum ra- 
dius of 6>^ ax = 0.8° and 0.7° for PSRs J0248+6021 and 
J2240+5832, respectively, reduces the background at low 
energies. We also remove events from nearby 1FGL sources 
within < 3° using the same bounded, energy- dependent ra- 
dius as applied to the pulsar. 

The on-pulse sp ectra were obtained with a maximum 
likelihood analysis (jMattox et alJll996j) of the LAT data 
within 15° of the pulsar, using the Fermi science tool "gt- 
like" . The likelihood method weights events from the target 
and background sources according to the energy-dependent 
PSF, which explains the need for such a large ROI. We ex- 
cluded time intervals when the ROI intersected the Earth's 
limb. In this paper, we modeled the pulsars using the func- 
tional form 



dN 
dE 



= N E 



exp 



cm^s^MeV -1 , 



(2) 



differential flux normalization 
T the photon index, and E c 



where Nq is the 
(phcm- 2 s- 1 MeV" 1 ), 
the cutoff energy. A pure power law (/? = 0) describes 
many gamma-ray sources, such as active galactic nuclei 
(lAbdo et al.ll2010al). T he 46 gamma-ray pulsars discussed 



Abdo et all (|2010bh are generally well-described by a 
simple exponential cutoff, j3 = 1, a shape predicted by 
outer magnetosphere emission models (see the Discussion, 
below). Models where gamma-ray emission occurs closer 
to the neutron star can have sharper "super-exponential" 
cutoffs, e.g. j3 = 2. 

The Galactic diffuse emission was modeled using the 
glLiem-v02 map cube based on six Galactocentric "ring" 
maps of N(H/) and W(CO) and on the spatial distribution 
of the inverse Comp ton intensity modeled by GALPROP 
(jStrong et al.l l2004 a.b). Diffuse extragalactic gamma-ray 



Gamma-ray data, analysis software, rotation ephemerides, 
and the diffuse background models are publically avail- 
able at the Fermi Scien ce Support Center, FSSC, 
http:/ / fermi. gsfc. nasa. g ov /ssc/ 1 



emission and residual instrument backgrounds were mod- 
eled jointly by the isotropic component isotropicJ,em_vQ2. 
Both models are available from the FSSC. All 1FGL sources 
within 20° were modeled with a power-law. Sources farther 
than 5° from the target pulsar were assigned fixed power 
law spectra, with parameters taken from the 1FGL source 
catalog. Spectral parameters for the pulsars and sources 
within 5° of the pulsar were allowed to vary. Systematic 
uncertainties were estimated by reapplying the fitting pro- 
cedures using bracketing IRFs where the effective area was 
shifted (linear extrapolations in log space) by ± 10% at 0.1 
GeV, ± 5% near 0.5 GeV, and ± 20% at 10 GeV. 

5.1. PSR J0248+6021 

PSR J0248+6 021 is called 1FG L J0248.3+6021 in the 
LAT catalog (|Abdo et al.l l2010al ). It is in the Galactic 
plane, 1.3° away from the bright gamma-ray source 
1FGL J02 40. 5+6113, identifie d as the Be star binary LSI 
+61° 303 (jAbdo et al.ll2009dl) . No other 1FGL sources are 
within 3°. The hig h E pulsar PSR J0 205+6449 in 3C 58 
is about 6° away (|Abdo et al.l l2009bh . PSR J0248+6021 
was discovered towards the end of EGRET's lifetime, 
and ephemerides contemporaneous with only a fraction of 
EGRET data exist. Before the Fermi mission, we searched 
for gamma-ray pulsations in all EGRET data for a range 
of {y, v) values extrapolated from our radio ephemeris. 
However, no pulsed signal was detected presumably because 
of the poor S/N of the EGRET data. 

Fermi LAT data for PSR J0248+6021 do however ex- 
hibit clear pulsations for the 3810 photons with energies 
> 0.1 GeV remaining after cuts. Fig. [1] shows the folded 
light curve of these events for different energy bands, along 
with the radio profiles used to derive the timing parameters. 
The light curve in each energy interval consists of one broad 
peak between 0.20 and 0.55 in phase, which dominates the 
highest energy band (> 1 GeV) but is progressively less pro- 
nounced with decreasing energy. The pulsation significances 
based on the bin- independent -ff-test (|de Jaeer fc Biischind 
120101 ) are ~ la, 6.5c, 8.2er, respectively, for the three en- 
ergy bands 0.1 - 0.3, 0.3 - 1.0, and >1 GeV. Between 0.3 
and 1 GeV, a structure with statistical significance of ^3 er 
(signal / ^background) appears between phases 0.6 and 0.7. 
However, it is not observed in the other energy bands, and 
a full likelihood analysis in that phase interval yields an ex- 
cess with a significance of only 2 a between 0. 1 and 100 GeV. 
We fit the main 7-ray peak (0.2 < <\> < 0.55) above 0.1 GeV 
with two half-Lorentzian functions. This shape provides a 
higher value of % 2 than others and accommodates differ- 
ences between the leading and trailing edges. The fit peak 
is at phase 6 = 0.39 ± 0.02 after the maximum of the ra- 
dio peak, with a full width at half maximum (FWHM) of 
0.20 ±0.02. The uncertainty in S is statistical. The bias due 
to the DM uncertainty in extrapolating the radio TOA to 
infinite frequency is negligible. The peak position does not 
vary with 7-ray energy within the statistical uncertainties, 
and the highest energy photon has 5 GeV at phase 0.31. 

We fit the on-pulse (0.2 < <j> < 0.55) data for 
PSR J0248+6021 with an exponentially cutoff power law 
(/3 = 1). The two neighbors mentione d above, as well as 
PSR J0218+4232 (jAbdo et al.l l2009alh were modeled in 
similar ways. Softer or sharper cutoffs (/? 7^ 1) yield essen- 
tially the same statistical significance. The likelihood ratio 
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test (2A(loglike), twice the difference of the logarithm of 
the likelihood maximum value) prefers the simple exponen- 
tially cutoff power law (j3 = 1) for PSR J0248+6021 to a 
simple power law (j3 = 0) by 8a. The spectral results for 
@ = 1 are listed in Table [TJ along with the integral photon 
flux Fioo and integral energy flux Gioo above 0.1 GeV. The 
first errors are statistical and the second are systematic. 
Fig. IT21 shows both the overall fit between 0.1 and 9 GeV 
(solid lines) with (3 = 1, and spectral points from likeli- 
hood fits in each individual energy band using a power-law 
spectrum. 

We performed a maximum likelihood analysis for 50 
days before and after the 2009 March glitch epoch. No 
change in the flux was observed. We searched the off-pulse 
data for a point source in the energy band 0.1-100 GeV 
at the radio pulsar position. No signal was observed. After 
scaling to the full pulse phase, we derived a 95% confidence 
level upper limit on the flux of 1 x 10~ 8 cm~ 2 s _1 . If a 
PWN is associated with the pulsar, deep observations will 
be required to see it . 

5.2. PSR J2240+5832 

PSR J2240+5832 is in the Galactic plane 0.6° away 
from the radio-qui et gamma-ray pu lsar PSR J2238+5903 
(5i4oo < 7/iJy, iRav et aD l2010f) . and 3° away from 
the Vela- like PSR J2229+611 4 associated with the 
"Boomerang" PWN and SNR (lAbdo et all 120093) . No 
other 1FGL source lies within 3°. The high background 
from nearby PSR J2238+5903 is the reason for the choice 
of small value of 6»™ ax = 0.7°. A total of 1208 photons 
with energies > 0.3 GeV were selected. Fig. [5] shows the 
gamma-ray pulse profiles in different energy bands, phase- 
aligned with the peak of the 1.4 GHz radio profile observed 
at Nancay. A single peak between phases 0.55 and 0.63 
appears above 0.3 GeV only. We tried different apertures 
to search for a signal below 0.3 GeV, but no excess was ob- 
served at the peak location. The highest //-test significance 
is 6.3 a and occurs above 1 GeV. Fitting above 0.3 GeV with 
a Lorentzian function places the peak at 5 = 0.58 ±0.01 af- 
ter the radio maximum, with an FWHM of 0.11 ±0.02. The 
peak position does not vary with 7-ray energy within statis- 
tical uncertainties, and the highest energy photon (20 GeV) 
lies at <fi = 0.72, just outside the peak phase range. 

We determined the on-pulse ((f) — 0.50 — 0.75) spec- 
trum using the maximum likelihood analysis within a 
15° aperture. Both nearby pulsars PSRs J2238+5903 and 
J2229±6114 as well as PSR J2240±5832 were modeled by 
a simple exponentially cutoff power law (/? = 1 in Eq. 2). 
Results are listed in Table [TJ The source is three times 
fainter than the other pulsar, leading to large uncertainties 
in the spectral parameters. The spectral index and cutoff 
are consistent with those seen for most young gamma-ray 
pulsars. We explored different models and the exponential 
cutoff (j3 = 1) is preferred to a power law by only 2a. 
Fig. IT21 shows both the overall fit between 0.1 and 100 GeV 
(solid lines), along with the spectral points from power-law 
likelihood fits to each individual energy band. 

As for PSR J0248+6021, we searched the off-pulse re- 
gion for a possible PWN at the pulsar position. No sig- 
nal was observed. Scaling to the entire phase range (a 
full rotation) yields a 95% confidence level upper limit 
on the flux of 2 x 10~ 7 cm~ 2 s~ 2 . We note however that 
the MILAGRO significance map for the region surround- 



ing the "Boomerang" PWN and PSR J2229±6114 extends 
intriguingly towards the Galactic plane, reaching 4 .7a at 
the position of PSR J2238±5903 (jAbdo et all 12009 fl ). The 
Tibet air sho wer array similarly see s a 2.5(7 excess at 
that location (|Amenomori et al.ll2010l ). Both MILAGRO's 
and Tibet's angular resolutions are poor, and TeV PWN's 
are notoriously offset from the pulsars driving them. PSR 
J2238+5903 is more energetic (£ = 9 x 10 35 erg s" 1 ) than 
PSR J2240+5832 and presumably closer, given that it is 
brighter (i*ioo = 6.8 ± 1.5 x 10~ 8 photons cm~ 2 s _1 above 
100 MeV). The TeV excess might be a conflation of two ad- 
jacent sources. Deeper observations with a Cherenkov im- 
ager array may be interesting. 

6. Discussion 

6.1. Gam ma- ray luminosity 

Armed with the gamma-ray integral energy flux Gioo and 
the distances d, we can now evaluate the lum inosity L 7 (Eq. 
1). We set /n = 1 as in lAbdo et al.l (|2010bD . and see below 
that the emission models in any case yield fn — 1. Defining 
G11 = 10 -11 erg cm~ 2 s _1 and d\ = 1 kpc = 3.1 x 10 21 cm 
gives L 7 = 1.2xl0 33 /o(Gioo/Gn)(d/di) 2 ergs" 1 . For PSR 
J0248+6021, summing the distance and Gioo uncertainties 
in quadrature yields L 1 = (1.4 ± 0.3) x 10 34 /si erg s — 1 . For 
PSR J2240±5832, the distance ambiguity leads to larger 
uncertainties: if the pulsar is in the outer arm, we obtain 
L 1 = (7.9 ± 5.2) x 10 34 in erg s"\ but if there is some 
unseen cloud along the line-of-sight creating an overdensity 
of electrons, and the pulsar is in the nearer Perseus arm, 
then L 1 = (2.2 ± 1.7) x 10 34 /^ erg s -1 . The corresponding 
efficiencies i] — L 1 /E are listed in Table 1. 

We compare these luminosities with the rule-of-tfmmb 
for young pulsars illustrated in Fig. 6 of the Fermi Pulsar 

Catalog (lAbdo et alJl2010bD . = 10 33 ^S/10 33 erg/s, 

where h stands for "heuristic" . For the E of our two pul- 
sars, L 7 is 1.5 x 10 34 erg s _1 , very near L 1 obtained 
for PSR J0248+6021. The measured value of L 1 for PSR 
J2240+5832 overlaps the rule-of-thumb for the distance of 
the Perseus arm, while the outer arm value is within the 
spread of values observed for the other pulsars. It is im- 
portant not to assign a distance to the pulsar based on a 
comparison of observed L 7 with expectations, since that 
may bias future attempts to use population modeling to 
refine emission models. 

6.2. Emission models 

These two pulsars have atypical gamma-ray pulse profiles : 
as for the six pulsars studied bv IWeltevrede et al.l (|2010t) . 
they have only one gamma-ray peak, and like four of those, 
they have a large offset from the radio pulse. Most gamma- 
ray pulsars have two peaks, the first lagging the radio beam 
by < 0.2 rotations. Hence, these two pulsars help extend 
the parameter space over which we can test gamma-ray 
emission models. 

To pin down the geometrical angles a and £, 
beyond the RVM work described above, we follow 
iRomani fc Wattersl (|2010f ) in modeling the gamma-ray 
pulse profiles. Byproducts are the flux correction factors 
/o, which reduce the uncertainty in L 7 . We test two sim- 
ple versions of outer magnetosphere pulse models. The first 
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model ('outer gap', or OG, 1 oman H ll99fih has emission 
starting at the 'null charge' surface and extending to the 
light cylinder. The second model ('two pole caustic', or 
TPC, iDvks fc Rudakl [2003) has emission starting at the 
star surface and extending to a perpendicular distance of 
0.75Rlc from the rotation axis or a radial distance of Rlc, 
whichever is less. In both cases, the emission comes from a 
zone spanning a characteristic fraction w — L^/E — 0.075 
of the open zone near the last-closed field lines. For PSR 
J0248+6021, we choose a linear intensity gradient across 
the gap but the modeled light curve shape is not very sen- 
sitive to the illumination across the gap zone, and narrow or 
uniformly illuminated gaps give similar results. In contrast, 
the very sharp 7-ray pulse of PSR J2240+5832 seems to 
require a relatively small range of field lines near the max- 
imum w to dominate the emission: here we used a simple 
Gaussian weighting, peaked around w — 0.075. We also plot 
the fits to the numerical plasma-filled (force free) model of 
iBai fc Spitkovskvl (|2010h . This so-called 'separatrix layer' 
(SL) model has a dense plasma and currents and posits 
emission from field lines extending well outside the light 
cylinder into the wind zone. 

Figs. and [5] show the goodness-of-fit surfaces in 
the (a, () plane for the two pulsars. "Goodness-of- 
fit" matching of the observed and modeled post-profiles 
uses the exponent i ally-t apered X3 weighting defined in 
IRomani fc Watterel pOlCh . a more robust test statistic than 
a simple % 2 . The regions selected by the RVM fits to the 
radio data are superimposed. In both cases, emission from 
near the light cylinder (as used in this OG model) seems 
necessary to provide good fits. In addition, emission from 
below the 'null charge' surface appears to create too many 
pulse components and too much off-pulse emission. For PSR 
J0248+6021, the best-fit curves in the radio-allowed region 
are near a = 46°, ( — 52°. This model produces a broad, 
merged peak at the correct phase with a good fit statistic 
X3 — 2.4. The inferred flux correction factor is fa = 1.06. 
The best fit TPC model in the allowed region is near the 
same angles but gives a substantially poorer fit, because 
of a second pulse component at 4> ~ 0.1. Here we have 
fa = 0.86. 

A similar exercise for PSR J2240+5832 shows rather 
similar regions of best fit for the two models. However, 
again the OG model provides a better fit (near a — 101°, 
C = 117°) with X 3 = 4.1. The poorer TPC fit at this po- 
sition has X3 = 7.5. In this case, the OG light curve is 
dominated by a single narrow pulse near <fi = 0.6. A much 
weaker component near <f> = 0.2 is faint or missing in the 
present gamma-ray light curve. Although these same two 
components are present in the TPC curve, the first peak 
after the radio pulse is comparably bright to the second, 
and for these a, £ angles significant emission is predicted 
for all pulsar phases. These differences indicate that, us- 
ing the parti cular realizations of the O G and TPC models 
described bv IRomani fc Wattersl (|2010h . the data are most 
consistent with the OG-type picture. For both fits, /n ~ 1. 
The SL models provide X3 At levels intermediate between 
the OG and TPC models. The relatively good fits tend to 
lie along a narrow band and the intersection with the RVM- 
allowed angles is similar to that proposed in the OG case. 

Thus, while both of these pulsars show single dominant 
components in the GeV light curves, the pulse profile mod- 
eling suggests a different origin. For PSR J0248+6021, the 
single broad component near phase 4> ~ 0.3 — 0.45 may be a 



merged double from a grazing sweep across the hollow 7-ray 
beam. In contrast, the very narrow pulse at phase <f> = 0.6 
for PSR J2240+5832 implies a cut through a caustic sur- 
face for a component normally identified with the second 
peak in young pulsars: the first peak is weak or absent. 




Energy(GeV) 




Energy(GeV) 



Fig. 12. On-pulse gamma-ray spectral energy distributions 
(SEDs) obtained with the Fermi Large Area Telescope. 
Top: PSR J0248+6021, in the phase range 0.2 < (f> < 0.55. 
Bottom: PSR J2240+5832 in the phase range 0.50 < <fi < 
0.75. Plotted points are from power-law likelihood fits to 
individual energy bands with > 2a detection significance 
above background for two degrees of freedom, otherwise an 
upper limit arrow is shown. The solid black line shows the 
maximum likelihood fit to a power law with exponential 
cutoff (Eq. 2). The dashed lines are ±ltr uncertainties in 
the fit parameters. 



7. Conclusions 

The Nangay survey of the northern Galactic plane led to the 
discovery of two remarkably similar young pulsars. They 
were subsequently timed for twelve years, during which pe- 
riod instrument precision improved steadily. The extensive 
timing data allowed proper motion determinations, and re- 
vealed a slowly varying DM for one pulsar. These data, to- 
gether with RM measurements obtained with polarization 
data, proved crucial in constraining the distance of PSR 
J0248+6021 to 2 kpc, much less than the value deduced 
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from the DM alone. The timing data also facilitated inves- 
tigations of the pulsar's environment. Radio polarization 
data constrained the pulsar geometry. 

Both pulsars have large spin-down power and are among 
the growing number of gamma-ray pulsars detected with 
the Fermi Large Area Telescope. Yet they are in the mi- 
nority of gamma-ray pulsars displaying only a single peak, 
with a large offset from the radio pulse. Gamma-ray emis- 
sion modeling constrains the regions where the electrons 
radiate at high energies. The gamma-ray luminosities ob- 
tained reinforce the growing evidence for how emission de- 
pends on spin-down power. 

We have studied radio, optical, X-ray, and TeV data for 
evidence of either accompanying pulsar wind nebulae or 
a bright optical pulsar. None was seen, but we presented 
evidence to encourage deeper X-ray and TeV searches. 
Sustained, accurate pulsar timing combined with multi- 
wavelength data are thus allowing progress in the under- 
standing of the emission mechanisms at work around neu- 
tron stars, and of the diffuse Galactic radiations at all en- 
ergies. 
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Table 1. Measured and derived parameters for pulsars PSR J0248+6021 and PSR J2240+5832. 



Pulsar name 


J0248+6021 


J2240±5832 


Right ascension, (J2000) 


02:48:18.617(1) 


22:40:42.939(4) 


T~\ 1* j_ ■ / ~f c\r\r\r\\ 

Decimation, (J 2000) 


+60:21:34.72(1) 


1 CO OO. O T f C\ ( A\ 

±58:32:37.52(4) 


Tl ' J- ' 1 /AT TT\\ 

Position epoch, (MJD) 


54000 


54000 


Galactic longitude, I (degrees) 


136.9 


106.6 


Galactic latitude, fe(degrees) 

Pulse frequency, v (s _1 ) 


0.698 


n 111 
-0.111 


4.6063 


7.1462 


first derivative oi pulse irequency, v (10 s ) 


— 1.169 


—0.7843 


Proper motion in right ascension, fi a (masyr -1 ) 


48 ± 10 


-6.1 ± 0.8 


Proper motion in declination, [is (masyr -1 ) 


48 ± 4 


-21.0 ± 0.4 


Timing data span (MJD) 


54682.7 to 55160.9 


54682.7 to 55202.6 


Rms timing residual (/is) 


181 


66 


Radio pulse width at 1.4 GHz, W50 (phase) 


0.12 


0.07 


Dispersion measure, DM (cm _3 pc) 


370 ± 1 


263.50 ± 0.05 


DM epoch, (MJD) 


54000 


54000 


Rotation measure, RM (rad m~ 2 ) 


-158 ±6 


24 ± 4 


Mean flux density at 1400 MHz, S1400 (mjy) 


13.7 ±2.7 


2.7 ± 0.7 


Radio spectral index 


-1.2 ± 0.1 


— 1.8 ± 0.3 


Radio-gamma-ray peak offset, 8 


0.39 ± 0.02 


U.00 ± U.U1 


Gamma-ray spectral Parameters 




Differential photon flux normalization, JVo (cm s 1 MeV 1 ) 


1.5 ±0.7 ±0.1 


n 1 _i_ n nfi _j_ a i a 
0.13 ± 0.06 ± 0.10 


Integral photon flux, F100 (> 100 MeV, 10~ 8 ph cm~ 2 s _1 ) . . 


4.0 ±0.4 ±0.5 


1 r . i ^ i n a 


Integral energy flux, G100 (> 100 MeV, 10~ n erg cm~ 2 s _1 ) . 


2.9 ±0.3 ±0.2 


1 A _l_ A A _L A O 

1.0 ± 0.4 ± 0.2 


Exponential cut-off energy, E c (GeV) 


1.2 ±0.2 ±0.1 


T _L A A _L 1 A 

5.7 ± 4.4 ± 1.0 


Power law index, T 


1.3 ±0.1 ±0.1 


1.8 ± 0.0 ±0.1 


Derived Quantities 






Spin-down power, E (erg s _1 ) 


2.13 x 10 35 


2 12x1 35 


Characteristic age, r c (kyr) 


63 


151 


Surface magnetic field strength, Bs (10 12 G) 


3.5 


1.45 


Magnetic field strength at the light cylinder, Blc (G) 


3150 


4860 


Distance from NE2001, d N E (kpc) 


> 43.5 


10.3i^ 3 


Distance used in this work, d (kpc) 


2.0 ±0.2 


7.7 ±0.7 or 3.8 ±0.8 


Luminosity, L 7 (E > 100 MeV, 10 34 erg s _1 ) 


1.4 ±0.3 


7.9 ±5.2 or 2.2 ± 1.7 


Efficiency, r) — L-y/E (%) 


7±2 


37 ± 25 or 10 ± 8 
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